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ABSTRACT 

We analyzed the longest phase-connected photometric dataset available for NGC 1313 
X-2, looking for the ~6 day modulation reported by Liu et al. (2009). The folded B 
band light curve shows a 6 day periodicity with a significance slightly larger than 3a. 
The low statistical significance of this modulation, along with the lack of detection in 
the V band, make its identification uncertain. 
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1 INTRODUCTION 

X-ray observations of nearby galaxies show a population of 
point-like, off-nuclear sources with luminosities (if isotropic) 
in excess of the classical Eddington limit for a IOMq com- 
pact object; they are referred to as Ultraluminous X -ray 
Sources (ULXs; see lFabbianoll2006l ; iFeng fe Soriall201ll for a 
review). Nowadays, hundreds of ULX candidates have been 
detected and many of them have been studied in detail 
(see e.g. iRoberts fc Warwick! l200d : IColbert fc Ptakl 12002 ' 



Swartz et all 12004 ILiu fc Bregmar] 120051; ILiu fc Mirabell 
20051 : IWalton et al.ll201ll : ISwartz et all 120111 ). Their high 



X-ray luminosities and short term variability suggest that 
the majority of these puzzling sources are likely to be 
accreting black holes (B Hs) in binary systems (see e.g. 
IZampieri fc Robertsll2009l and references therein). Although 
X-ray data alone have provided evidence that ULXs are dif- 
ferent from stellar-mass Galactic BHs - either a different 
class of BHs or a different accretion state - , X-ray spectral 
and timing results remain consistent with various alternative 
scenarios, characterized by BHs of different mass and ori- 
gin. These go from the challenging intermediat e-mass BHs 
of sa 10 2 - 1O 4 M (|Cofbert fc Mushotzkvlll999l '), to massive 
stellar BHs of ~ 25 — 80M© formed from the direct col- 
lapse of low-metallicity m assive stars (Ma pelli et al.l [20091 : 
IZampieri fc Roberts! [20091 ') . to stella r-mass BHs (< 20M G ) 
accre ting above the Eddington limit l|King et al.ll200l1 ; lKing| 
l200Sl h 

A necessary pre-condition for determining (or at least 
constraining) the BH mass is to measure the period of a 
ULX binary system. This represents a crucial preliminary 
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step, required to perform an efficient spectroscopic follow- 
up in search of radial velocity variations and hence mea- 
sure the mass function of the system, and so the BH mass. 
The best way to quantify the binary period is to study 
the optical emission of the counterpart and to determine 
its periodic modulation. The same strategy was success- 
fully adopted for measuring the orbital period and mass 
of the compact object in Galact ic BH X-ray binaries (e.g. 
Ivan Paradiis fc McChntocklll995h . 

To date a single dedicated monitoring campaign was 
performed on NGC 1313 X-2 with HST, providing the only 
tentative optical periodicity available. NGC 1313 X-2 is lo- 
cated in the outskirts of the bar red spiral galaxy NGC 1313 
at a distance of 3.7-4 .27 Mpc l|Tullv Ill98sl ; iMendez et~ail 
120021 : iRizzi et al.ll2007l ) . Its observed X-ray luminosity varies 
between a few xlO 39 erg/s and ~ 10 40 erg/s in the 0.3-1 
keV band (e.g. lMucciarelli et al.ll2007l ; lFeng fc Kaaretll2006t ). 
The source has bee n extensively studied in the X-ray and 
optical bands (e.g. |Zampieri et al. 1 12004 IMucciarelli et al.l 
120071 : iGrise et al.l 20081 and references therein). The large 
amount of data available make this object a cornerstone for 
the study of ULXs. It belongs to a handful of UL Xs clearly 
associated with stellar optical counterparts (e .g. ILiu et all 
|2004 IMucciarelli et al.ll2005l : ISoria et al.ll2005l h These opti- 
cal sources appear to be al most ubiquitously hosted in young 
stellar environments (e.g. iPakull et al.ll200a : I Ramsey et al.l 
120061 : ILiu et all l2007f) and have properties consistent with 
those of young, massive stars. However, some ULXs appear 
to be associated with older stellar populations and at least 
one possible later type stellar counterpart is now k nown 
jFeng fc KaaretJ 120081 : IRoberts, Levan fc Goadl [2008), al- 
though its spectral classification may be affected by sig- 
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nificant galactic and extra-galactic reddening (I Grise et al.l 
l2006h . 

The optical counterpart of NGC 1313 X-2 was first 
identified on an ESO 3.6 m R band image thanks to 
Chandra's accurate astrometry, after the X -ray image was 
regist ered on the position of SN 1978K (lZampieri et al.l 
l2004h . ESO Very Large Telescope (VLT) images of the 
field showed that the counterpart was actually composed 
of two distinct objects (CI and C2), separated by ~ 0.7" 
l|Mucciarelli et al1l2005l ). Further refinement in the astrom- 
etry and accurate modelling of the optical emission indicated 
that o bject CI was the more l ikely counterpart (Liu et"al] 
120071 : IMucciarelli et all l2007t IPatruno fc Zampieril EToOSl) . 
However, this was established beyond any doubt by the 
detection of the HeIIA4686 emission line in its optical 
spectrum, a chara c teristic imprint o f X-ray irradiation 
jPakull et all 120061 : iGrise et all 120081 ). This star has an 
extinction-corrected absolute magnitude Mb ~ —4.5 mag 
and c olors (B - V) ~ -0.15 mag and (V - I ) ~ -0.16 
mag l|Mucciarelli et al.l [20071 : IGrise et al.l 12008) 1 . consistent 
with a B spectral type. 

The stellar environment of NGC 1313 X-2 has also 
provided interesting constraints. There are two groups of 
young stars spread out over ~ 200 pc. Isochrone fit- 
ting of the colour-magnitude diagram of these groups has 
been attempted an d provides clus t er ag e s of 20±5 Myrs 
(iPakull et all 120061 ; iRamsev et ail l200d : iLiu et all 120071 ; 
IGrise et alj|2008| ). As several other ULXs, NGC 1313 X-2 is 
also associated with a very extended (~400 pc) optical emis- 
sion nebula that gives import ant information on the energet- 
ics and lifetime of the system (|Pakull et al.ll2002T ) . Assuming 
that it is formed in one or more explosive events or that its 
mechanical energy comes from the ULX wind /jet activity, 
the characteristic age and energetics of the nebula turn out 
to be ~1 Myr and w iq 52 ^ 53 e rg or ~ 4 x 10 39 erg/s, re- 
spectively ((Pakull et al.| [2006). The most up-to-date binary 
model calculation, including X-ray irradiation effects, finds 
consistency between all the available optical measurements 
if CI is a terminal-age main sequence or early giant donor of 
12 — 15M0 undergoing Roche-lobe overflow. The same cal- 
culations provide estimates of the ma sses also for the BH, 
which is in bet ween 20 and ~ lOOMm (IPatruno fc Zampieril 
|2010| ; see also ICopperwheat et all |2007l ) for a similar re- 
sult including star+disc irradiation but without considering 
binary evolution effects). 

Quite recentlv. lLiu et all (|2009T ) found a possible period- 
icity of 6.12±0.16 days in a monitoring campaign of the opti- 
cal counterpart of NGC 1313 X-2 performed with the Hubble 
Space Telescope (HST). This modulation was interpreted as 
the orbital period of the binary system. Three cycles were 
detected in the B band, while no modulation was found in 
V . Previous studies carried out on the available HST an d 
VLT observations led to negative results (|Grise et al] |2008). 
More recently, lack of significant photometric variability on 
a new sequ e nce o f VLT observations has been reported by 
IGrise et all (|2009l ). In principle, the detection of the orbital 
period would definitely confirm the identification of the op- 
tical counterpart and the binary nature of this system. Most 
importantly, it would open the way to perform a dynamical 
measurement of the BH mass. 

Here we present a re-analysis of the joint VLT+FORS1 
and #ST+WFPC2 photometric observations of NGC 1313 



X-2 obtained during the years 2007-2008, with the aim of 
clarifying the st atistical s i gnific ance of the orbital periodic- 
ity identified bv lLiu et al.l (|2009l ). We did not consider previ- 
ous VLT observations taken in 2003-2004 because they can- 
not be phase-connected to the more recent ones. In §[2] we 
present the data reduction procedure that we have adopted, 
while in §[3] we show the results of the statistical analysis. 
§ [4] summarizes our conclusions. 



2 VLT AND HST OBSERVATIONS 

NGC 1313 X-2 was observed with VLT +FORS1 between 
October 2007 and March 2008 (11 epochs; [Crise et alj|200_h 
and wit h j.5T+WFPC 2 between May and June 2008 (20 
epochs ; ILiu et all 120091 ) . During the VLT observations the 
sky was clear and the average seeing was in the range 
0.8"-1.2". The quality of the WFPC2 images is fair, despite 
the degradation of the central PC chip. A log of the obser- 
vations is reported in Table [1] We re-analyzed the whole 
dataset in a homog eneous way, look ing for the ~6 day peri- 
odicity reported bv lLiu et al. [ ||2009h . 

Two exposures were taken a few minutes apart each 
night. For the HST dataset, we used the calibrated data 
from the WFPC2 static archive. After performing standard 
image reduction in the IRAF environment, the exposures 
were combined together and cleaned for cosmic rays. To ac- 
curately photometer the objects, we use d AIDA (Astronom- 
ical I mage Decomposition and Analysis; lUslenghi fc Falomol 
2008), an IDL-based package originally designed to per- 
form two-dimensional point-spread-function model fitting of 
quasar images. For the analysis of the WFPC2 exposures, 
we loaded into AIDA the appropriate point-spread-function 
simulated with Tiny Tim v. 6.33- As the background did 
not vary significantly among the different HST exposures, 
we decided to keep it fixed at the average value computed 
from all the observations. 

Analyzing VLT and HST measurements together re- 
quires attention to be paid to the systematic differences 
between the two photometric systems. We then first con- 
verted the HST instrumental magnitudes to the standard 
UBVRI photometric system using the updated transforma- 
tion equations and coefficients published in Dolphin (2009) 
for the appropriate instrumental gain (which is equal to 7 in 
our casejj. The color correction te rm was computed adopt- 
ing the (B — V) color reported in IMucciarelli et al.l l|2005l) . 
In spite of this, residual systematic differences between the 
two photometric systems might still be present and affect 
our measurements. In particular the color correction term 
is sensitive to the overall bandpass (telescope plus atmo- 
spheric response) of the instruments. We then decided to 
perform differential photomet ry of the target with r espect to 
a nearby field star (star D in lZampieri et al1l2004l ). located 
on the same chip in both instruments. The coordinates and 
average BVR magnitudes of the optical counterpart of NGC 
1313 X-2 (o bject CI) and the refere nce star (object D) are 
reported in IMucciarelli et al.l l|2005l ). The reference star is 
brighter than the target and has a low root mean square 
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Table 1. Log of the VLT+FORS1 and ffi?T+WFPC2 photometric observations of the field around NGC 1313 X-2, along with the B 
and V band differential photometry of object CI with respect to a reference field star (see text for details). 



Obs. Date MJD Exposure Telescope+Instr. AB a B AV a v 



00 


1 


2007-10-21 


54394 


.362587 


242x2 


VLT+FORS1 


3. 


.140 


0. 


069 








2 


2007-11-15 


54419 


.220472 


242x2 


VLT+FORS1 


3. 


.220 


0. 


045 








3 


2007-11-15 


54419 


.277624 


242x2 


VLT+FORS1 


3 


.261 


0. 


076 








4 


2007-11-16 


54420 


.255469 


242x2 


VLT+FORS1 


3 


.220 


0. 


049 








5 


2007-12-06 


54440 


.075805 


242x2 


VLT+FORS1 


3. 


.172 


0. 


042 








6 


2007-12-10 


54444. 


.173203 


242x2 


VLT+FORS1 


3. 


169 


0. 


.052 








7 


2007-12-14 


54448 


.121299 


242x2 


VLT+FORS1 


3 


.236 


0. 


071 








8 


2008-01-31 


54496 


.073081 


242x2 


VLT+FORS1 


3. 


.162 


0. 


058 








9 


2008-03-02 


54527. 


.059338 


242x2 


VLT+FORS1 


3. 


.215 


0. 


066 








10 


2008-03-05 


54530. 


.053776 


242x2 


VLT+FORS1 


3. 


.172 


0. 


.053 








11 


2008-03-08 


54533. 


.056052 


242x2 


VLT+FORS1 


3 


.244 


0. 


065 








12 


2008-05-21 


54607. 


.911122 


500x2 


HST+WFPC2 


3 


.217 


0. 


.075 










2008-05-21 


54607. 


.927094 


400+700 


FST+WFPC2 










4.759 





.043 


13 


2008-05-22 


54608. 


.043761 


500x2 


FST+WFPC2 


3. 


.098 


0. 


.072 










2008-05-22 


54608. 


.059733 


400+700 


FST+WFPC2 










4.824 


0. 


.045 


14 


2008-05-23 


54609. 


.042372 


500x2 


HST+WFPC2 


3 


.058 


0. 


071 










2008-05-23 


54609. 


.058344 


400+700 


FST+WFPC2 










4.816 





.044 


15 


2008-05-24 


54610. 


.040983 


500x2 


F5T+WFPC2 


3 


.192 


0. 


.074 










2008-05-24 


54610. 


.056955 


400+700 


HST+WFPC2 










4.841 


0. 


.045 


IB 


2008-05-25 


54611. 


.038900 


500x2 


HST+WFPC2 


3 


.190 


0. 


.074 










2008-05-25 


54611. 


.054872 


400+700 


FST+WFPC2 










4.855 





.046 


17 


2008-05-26 


54612. 


.104178 


500x2 


FST+WFPC2 


3. 


.241 


0. 


.075 










2008-05-26 


54612. 


.120150 


400+700 


F5T+WFPC2 










4.917 


0. 


.047 


18 


2008-05-27 


54613. 


.179178 


500x2 


F5T+WFPC2 


3 


.226 


0. 


.075 










2008-05-27 


54613. 


.245150 


400+700 


FST+WFPC2 










4.815 





.045 


19 


2008-05-28 


54614. 


.178483 


500x2 


FST+WFPC2 


3. 


.185 


0. 


.074 










2008-05-28 


54614. 


.243761 


400+700 


HST+WFPC2 










4.877 


0. 


.046 


20 


2008-05-29 


54615. 


.177789 


500x2 


H5T+WFPC2 


3 


.063 


0. 


.072 










2008-05-29 


54615. 


.243066 


400+700 


FST+WFPC2 










4.850 





.045 


21 


2008-05-30 


54616. 


.097233 


500x2 


FST+WFPC2 


3. 


.110 


0. 


071 










2008-05-30 


54616. 


.113205 


400+700 


FST+WFPC2 










4.775 


0. 


.043 


22 


2008-05-31 


54617. 


.109039 


500x2 


HST+WFPC2 


3 


.178 


0. 


.075 










2008-05-31 


54617. 


.170150 


400+700 


FST+WFPC2 










4.854 





.046 


23 


2008-06-01 


54618. 


.692372 


500x2 


FST+WFPC2 


3. 


.326 


0. 


077 










2008-06-01 


54618. 


.708344 


400+700 


FST+WFPC2 










4.762 


0. 


.043 


24 


2008-06-02 


54619. 


.173622 


500x2 


HST+WFPC2 


3 


.200 


0. 


.075 










2008-06-02 


54619. 


.238205 


400+700 


FST+WFPC2 










4.835 





.045 


25 


2008-06-03 


54620. 


.105567 


500x2 


FST+WFPC2 


3. 


.189 


0. 


.074 










2008-06-03 


54620. 


.166678 


400+700 


HST+WFPC2 










4.788 


0. 


.043 


26 


2008-06-04 


54621. 


.171539 


500x2 


HST+WFPC2 


3 


.071 


0. 


.072 










2008-06-04 


54621. 


.236122 


400+700 


FST+WFPC2 










4.899 


0. 


.046 


27 


2008-06-05 


54622. 


.104178 


500x2 


FST+WFPC2 


3. 


.101 


0. 


.072 










2008-06-05 


54622. 


.164594 


400+700 


HST+WFPC2 










4.888 


0. 


.046 


28 


2008-06-06 


54623. 


.102789 


500x2 


HST+WFPC2 


3 


.316 


0. 


077 










2008-06-06 


54623. 


.163206 


400+700 


FST+WFPC2 










4.913 





.046 


29 


2008-06-07 


54624. 


.034733 


500x2 


FST+WFPC2 


3. 


.419 


0. 


080 










2008-06-07 


54624. 


.091678 


400+700 


FST+WFPC2 










4.872 





.045 


30 


2008-06-08 


54625. 


.100706 


500x2 


FST+WFPC2 


3 


.257 


0. 


.075 










2008-06-08 


54625. 


.161122 


400+700 


HST+WFPC2 










4.899 


0. 


.047 


31 


2008-06-09 


54626. 


.611817 


500x2 


FST+WFPC2 


3. 


.285 


0. 


.075 










2008-06-09 


54626. 


.627789 


400+700 


ffi>T+WFPC2 










4.833 





.045 



variability (~ 0.05 mag in the VLT and < 0.02 mag in the 
HST exposures). For similar reasons, during all the obser- 
vations performed with /f5T+WFPC2, the field was always 
oriented in the same direction and the target and reference 
star were always located on the same position on the central 
PC chip. 

Figures [1] and [2] show the B and V band light curves of 
NGC 1313 X-2 obtained in this way. The measured values of 



the differential magnitudes (AB = Bu — Bci, AV = Vu — 
Vci) and the corresponding errors (o"s, av) are reported in 
TableQ] The mean values and standard deviation of different 
samples of data are reported in Tabled A comparison of our 
B and V band light curves with those of iLiu et al.l (120091 ) 
(shifted by a constant value) shows that they are consistent 
within the errors, apart from the V band measurements of 
two observations (n. 24 and 31) that differ by 0.1 mag. 



4 L. Zampieri, D. Impiombato, R. Falomo, F. Grise and R. Soria 




105110 115120 125130 135 140 145 
MJD-54390 



220 



225 230 
MJD-54390 



235 



Figure 1. Joint VLT+HST light curve of NGC 1313 X-2 in the B band. The squares are the VLT+FORS1 data, while the circles 
represent the HST+ WFPC2 observations. The magnitudes are the difference with those of a reference field star. The data cover a period 
of ~ 7.5 months. The solid {cyan) line is the best fitting sinusoid with a period P = 5.68 days. 

Table 2. Average magnitudes, standard deviation and best fitting parameters of sinusoidal fits to different samples of data 



Sample 


< AS > 


< AV > cr* 


pa 


A a 




X 2 (dof) 








(days) 








~VLT+HST b (B band) 


3.198 


0.078 


5 68+ 01 
°- D8 -0.0l 


06+ 02 


27 0+29° 
Zl -31° 


21.0 (28) 








6 01+ 01 
°- ul -o.oi 


06+ 03 


q^O+25 
9t) _ 23 o 


21.9 (28) 


VLT (B band) 


3.201 


0.038 






86°+^ 




HST (B band) 


3.196 


0.094 


6 17+ ' 18 

-0.17 


10+ 002 

u - lu -0.02 


13.4(17) 



VLT+HST {V band) 



4.843 



0.047 



* Standard deviation 
a Errors are 90% confidence intervals. 
b Parameters of the two deepest minima. 



The same result is obtained by comparing the light curve of 
iLiu et all 1120091) with ours, before performing the UBVRI 
magnitude system conversion. 

Our data show clear short term (~ 1 day) variability. 
As can be seen from Figure Q3 the B band HST data have 
larger root mean square variability than the B band VLT 
ones. Similarly, comparing Figures [T] and [5] it appears that 
the B band light curve has larger variations than the V band 
data. 

Superimposed on these short term stochastic changes, 
the HST V band light curve does not show significant regu- 
lar variability (Figure [2]), whereas the HST B band dataset 



shows an approximately sinusoidal modulation with a period 
of 6 days (Figure [T]). The maximum peak-to-peak variation 
in the B-band light curve is 0.36 mag. 



3 RESULTS 



Following ILiu et al.l |2009l ). we fitted the HST and 
VLT+HST B-band light curves with a sinusoid: 



AB =< AB > +^sin(27r(t - ti)/P + . 



(1) 
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-4.70r 




-5.00 1 



220 



225 230 

MJD-54390 



235 



Figure 2. HST light curve of NGC 1313 X-2 in the V band. The 
magnitudes are the difference with those of a reference field star. 
The solid (cyan) line is a sinusoid with the same period and phase 
obtained from the fit of the joint VLT+HST B band data, but 
with a fixed amplitude of 0.04 mag. 



where A, P and <f) are the amplitude, period and phase, re- 
spectively, and ti — 54390 is a reference epoch. The best 
fitting parameters of the fit are reported in Table [2] Al- 
though the values of the period and phase of the HST and 
VLT+HST fits are not consistent at the 90% confidence 
level, they are in agreement at the 3-a level. Furthermore, 
the x 2 surface for the VLT+HST dataset has several (at 
least 5) pronounced minima in the range ~5. 7-7.1 days (5.7, 
6.0, 6.2, 6.7, 7.1 days, respectively), with similar values of the 
X 2 - The values reported in Table [2] refer to the two deepest 
minima. Considering these caveats, the fits appear to return 
consistent results. The value of the amplitude and period of 
the HST fit are in agreemen t, within the errors, with those 
reported bv lLiu et all i|2009r ). 

The error on the period for the VLT+HST light curve 
is such that AP(T/P) < 6 days, where T ~ 200 days is the 
interval between the first and last observation (sampling in- 
terval). This is not larger than P, indicating that the two 
datasets can be phase connected. Indeed, this represents the 
longest phase-connected photometric dataset available for 
NGC 1313 X-2. A fit of the VLT data alone was also at- 
tempted and returned P = 5.7 days, A = 0.04 mag and 
(f> — 28°. Computing the errors for one interesting parame- 
ter (while holding the others fixed) , all these values turn out 
to be consistent with those reported in Table[2] However, the 
amplitude of the sinusoidal fit (for fixed period and phase) is 
consistent with zero at the 90% confidence level. Therefore, 
it is not possible to infer evidence of periodicity from the 
VLT data alone, but they appear to be consistent with the 
sinusoidal modulation observed in the HST observations. 

We checked possible systematics effects induced by the 
transformation from the VEGAMAG to the UBVRI sys- 
tem fitting the function in equation (fjrj directly to the HST 
F450W band light curve. We found values for the period 
and amplitude completely consistent with those of the HST 
B-band light curve reported in Table[2](P = 6.21+°;^ days, 
A = O.llJ^a mag). The % 2 is 31.5 for 17 d.o.f., correspond- 
ing to a null hypothesis probability of 0.0174. The rather 
large value of the x 2 mav indicate that, at this level of 



accuracy, we may be sensitive to small deviations of the 
light curve from a perfect sinusoidal shape. The phase is 
4> = 169° ± 14°. This is clearly different from the phase re- 
ported in Table[2]because the reference epoch (MJD=54390) 
is ~ 215 days before the HST observations and small dif- 
ference in the period estimate (0.04 days) can accumulate 
to give rise to a phase shift of ~0.04 daysx (220/6) ~ 1.4 
days~ 80°. 

We tried to assess the statistical significance of the ap- 
parent B band modulation in two different ways. First we 
performed a Lomb-Scargle periodogram analysis of all the 
unevenly sampled time series, including both the VLT and 
the HST observations. The maximum power is at frequency 
of 1.93 x 10 -6 Hz, corresponding to a period of 6 days. How- 
ever, the null hypothesis probability is quite large (16.5%), 
meaning that the statistical significance of the modulation is 
low. In order to increase the signal-to-noise, we then decided 
to bin the light curve in M = 6 bin intervals and perform 
an epoch folding period search. The peak of the distribution 
is x 2 ~ LI .2 at P — 6 days. For v = M — 1 = 5 degrees of 
freedom, this corresponds to a significance level of 2.8er. No 
other significant peaks were found in the interval between 3 
and 9 days. As the validity of the y 2 dis tribution is limited 
to large samples, following iDaviesi ( 1990t ) we estimated the 
significance of the 6 day period using also the L statistics, 
which is statistically sound for all sample sizes. The peak 
of the distribution is L = 6.75 for v\ = M — 1 = 5 and 
V2 = N — M = 25 degrees of freedom, where N — 31 is the 
number of observations. This corresponds to a null hypoth- 
esis probability of 0.0004 or a significance level of 3.6<r. In 
Figure [3] we show the folded light curve along with its best 
fitting sinusoid (P = 6.0 days, A = 0.08, <f> = 54°). The 
parameters of the sinusoid are consistent with those of the 
second best fit of the unbinned light curve (see Table [5]). 

In contrast with the B band, the V band light curve 
shows a rather stochastic variability (although the first 6- 
7 observations appear to follow a sinusoidal behaviour). We 
tried a sinusoidal fit also of this dataset and found an accept- 
able minimum with a periodicity P = lO.Str' 2 , days that, 
however, is not statistically meaningful (significance from a 
Lomb-Scargle periodogram analysis ~ 56%). As shown in 
Figure [2] a sinusoidal modulation with the same period and 
phase obtained from the fit of the B band data and an am- 
plitude < 0.04 mag may be marginally consistent, within 
the errors, with the V band data (reduced x 2 1-4). 

Repeati ng a similar anal ysis on the F450W band mea- 
surements of lLiu et al.l |2009l) (their Table 1), we found that 
the maximum power in the Lomb-Scargle periodogram cor- 
responds to a period P — 6.13 days and to a null hypothesis 
probability of 16%. The peak of the x 2 distribution after 
binning and folding the light curve as described above is 
X 2 = 15.8 at P = 6.3 days. Adopting the L statistics, this 
corresponds to L = 13.8 for V\ = 5 and V2 = 14 degrees 
of freedom, or to a significance level of 4a (null hypothesis 
probability of 0.000054). This value is very close to that es- 
timated from our measurements, indicating that the main 
re ason for th e differ ence with the significance level reported 
bv lLiu et~a l, (2009) (6a) relies on the statistical treatment 
of the data and not on the different photometric analyses. 
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Figure 3. Binned light curve (6 bins) of the B band VLT+HST 
dataset of NGC 1313 X-2, folded over the best estimate of the 
period (6 days). The best fitting sinusoid is also shown. The phase 
is measured from MJD=54390. 

4 DISCUSSION AND CONCLUSIONS 

We re-analyzed the longest phase-connected photometric 
dataset available for NGC 1313 X-2 using VLT+FORS1 and 
BST+WFPC2 observations taken between October 2007 
and June 2008. B-band differential photometry with respect 
to a nearby isolated and non-s aturated field star confirms the 
6 day modulation detected bv lLiu et all (|2009l ) in the HST- 
F450W data alone. No significant periodic variability was 
found in the V band. Binning the B band light curve, the 
statistical significance of the 6 day modulation turns out to 
be slightly larger than 3cr. No other significant oscillations 
were found. 

Although much of the work presented in this paper 
is a re-reduction and analysis of previously published data 
(|Liu et al.ll2009l ; iGrise et alj|2009l ). we emphasize that there 
are two important differences in the approach that we fol- 
lowed here. We adopted differential photometry with respect 
to a reference field star (instead of absolute photometry) 
and, most importantly, we performed a full statistical anal- 
ysis of the significance of the modulation. Differential pho- 
tometry allows us to minimize the effects of absolute cali- 
bration uncertainties and leads to small differences in the 
flux measurements. All but two F-band measurements are 
in agreement with what reported previously in the litera- 
ture. However, the major improvement of our work consists 
in a careful statistical re-analysis of the data, that includes 
binning the li ght curve and using the L statistics for small 
sample sizes (|Davieslll990l ). While excluding the VLT ob- 
servations from the analysis does not affect significantly our 
conclusions, we find that they can be phase-connected to the 
HST dataset and are consistent with it. They also show that 
the optical luminosity of the counterpart is quite steady on 
timescales of months-years and that the accretion disc does 
not have a large stochastic variability, which is encouraging 
for photometric studies of this nature. At the same time, 
they also show that detecting periodic optical variability in 
ULX counterparts is possible both from space and from the 
ground, if 8+ meter class telescope are available. Necessary 
condition in the latter case is having optimal seeing and 
rather isolated counterparts. 



The observed optical emission of NGC 1313 X-2 origi- 
nates from the intrinsic emission and X-ray heating of the 
dono r star and the accretion disc (e.g. iPatruno fc Zampieril 
2008). From the maximum amplitude of the sinusoid that 
is marginally consistent with the VLT data, we infer a 
maximum rms sinusoidal variability of 0.03 mag for the 
VLT data, much smaller than that of the HST observations 
(~ 0.07 mag), indicating that there are different levels of 
optical activity whose origin is unclear. This is confirmed by 
the analysis of previous observations of NGC 1313 X-2 per- 
formed in 2003-2004, when the B band light curve showed an 
intermediate level of variability (standard deviation ~0.056 
mag) with respect to the two datasets considered here. There 
may be different reasons for this behavior. Some B stars have 
intrinsic variability on a timescale of a few hours that may be 
different from one cycle to the ne xt, even though the base- 
line is constant (e.g. p Cep stars; [Saesen et al"1 T2010). This 
variability may be superimposed to a short term stochastic 
variability from the disc, and appear random because the 
timescale is too short compared to the time resolution of 
our observations. However, if we consider the standard de- 
viation of the residuals with respect to the sinusoidal fit, 
the behaviour of the VLT and HST datasets is more sim- 
ilar (standard deviation ~0.05 mag for VLT, ~0.06 mag 
for HST). Therefore, it is likely that the extra variability 
observed in the B band HST data comes from the regular 
sinusoidal modulation. If it is so and if the changing view 
of the X-ray irradiated/non- irradiated donor surface gives a 
significant contribution to the modulation of the optical flux, 
an increased amplitude of the modulation would imply an in- 
crement in the X-ray heating and, in turn, in the average op- 
tical luminosity. We know that t he X-ray flux of NGC 1 313 
X-2 is indeed quite variable (e.g. iMucciarelli et al .1120071 ). A 
recent Swift monitoring campaign of the source, performed 
in 2009, seems to show a bimodal distribution of the fluxes 
which may in fact suggest rapid and recurrent X-ray flares 
(with variations of a factor of 3-4; Grise et al., in prepara- 
tion) , although the correlation with optical variability is not 
yet firmly established. 

At first glance, interpreting the extra HST variability 
as caused by increased X-ray heating does not appear to be 
consistent with the observed average AB reported in Ta- 
ble [5] which is not significantly different between the VLT 
and HST observations. On the other hand, as already men- 
tioned, it is important to remind that there may be possible 
residual systematic uncertainties in the conversion between 
the VEGAMAG and UBVRI photometric systems. In par- 
ticular, in performing the conversion, we applied a color cor- 
rection (to bot h objects CI and D) est imated from the col- 
ors reported in iMucciarelli et al.l |200S). However, we know 
that the color of object CI is variable and is likely to become 
bluer if X-ray heating increases. Furthermore, we found that 
the inferred (B — V) color of the reference object D in the 
HST data turns out to be sli ghtly redder (1.6 instea d of 
1.4 mag) than that reported in IMucciarelli et al.l l|2005l ) (al- 
though they are in agreement within the errors). Both effects 
would tend to change the average AB of the HST data and 
shift them upwards in Figure [T] as expected if X-ray irradia- 
tion has increased. For an error and/or variation in the colors 
of 0.15 mag, the average AB would vary by 0.035 mag, which 
is significant in comparison with the amplitude of the modu- 
lation. Even assuming that the average luminosity does not 
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vary much between the two datasets, one may still think 
to some particular conditions under which the modulated 
fraction of the emission varies without significant changes 
in < AB >. For example, when X-ray irradiation is higher, 
the disc may be partly obscured by material blown off from 
the inner regions. Or a change in the height of the outer ac- 
cretion disc may partly induce variations in the obscuration 
of the companion. While the total solid angle of disc+star 
seen by the source (and hence the fraction of intercepted X- 
ray photons) remains constant, the relative contribution of 
star and disc may change, inducing variations in the modu- 
lated fraction of the emission. Clearly no definite conclusion 
can be reached without further accurate observations. 

Also the smaller variability and absence of a detectable 
modulation in the V band is puzzling. Model calculations 
of the irradiated plus donor disc emission show that at 
longer optical wavelengths the donor spectrum declines more 
rapidly than the irradiated disc spectrum. Therefore, the 
contamination from the disc is comparatively stronger in 
the R and V band with respect to the B and U bands 
l)Patruno fc Zampieril I2OI0I ). As the emission from the ac- 
cretion disc is not orbitally modulated, any variation in the 
donor star emission would induce a slightly smaller change 
in the R and V bands than in the B and U bands. Perhaps 
the irregular variability comes more from the disc, while the 
phase-modulated variability more from the irradiated star. 
However, our analysis indicates an upper limit of 0.04 mag 
on the V band modulation, significantly lower than the in- 
ferred amplitude of the modulation in the B band (0.09 
mag). Furthermore, ellipsoidal modulatio ns are also likely 
to be important in these conditions (e.g. iBochkarev et al.1 
1 19791 ). with two maxima/minima per orbital revolution. If 
they dominate, no significant wavelength dependence of the 
modulation is expected and the observed 6 day periodic- 
ity would correspond to an orbital period of 12 days. In 
case both X-ray irradiation and ellipsoidal modulation con- 
tribute to the observed variability, the combined effect may 
be an asymmetric light curve with two minima of different 
depth (not easily detectable with the accuracy of present 
measurements). However, even in this case, the light curve 
may show a single maximum/minimum at superior/inferior 
conjunction, if the X-ray flux at the stellar surface is suffi- 
ciently high. Further observations and a detailed joint mod- 
elling of the irradiated accretion disc and ellipsoidally dis- 
torted donor are required to assess this point. We note that 
any signature of X-ray irradiation/reprocessing in the opti- 
cal data would provide evidence against beaming, as in that 
case the fraction of X-ray photons intercepted by disc and 
the star would be small. 

Another possibility for explaining the smaller V band 
variability may be the contribution of emission lines. The B 
band contains both the HeIIA4686 and H/3 lines, while the V 
band has only the continuum. If irradiation and/or orbital 
effects are causing the lines to vary in addition to the contin- 
uum, perhaps they may induce an extra B band variability. 
But the observed strength of the excess B band variabil- 
ity (~5% over the continuum) requires signifi cant line emis- 
sion a nd hence a very high equivalent width. iRoberts et al.l 
(2010) found that the equivalent width of the Hell A4686 
line displays variations between 2 Aand 11 A, not large 
enough to explain the observed variability. Finally, changes 
in the absorbing column towards the source may also con- 



tribute to induce differential variations in the B and V bands 
light curves (e.g. E(B — V) variability of ~0.01 mag would 
cause AB and AV variability of ~0.04 mag and ~0.03 mag, 
respectively). 

If the modulation is real and represents the orbital 
period of the ULX, it would place rather compelling con- 
straints on the pro perties of the system and als o on the donor 
and BH masses. iPatruno fc Zampieril (|2010h have shown 
that, imposing the 6 day periodicity and using all the other 
constraints available for this source from X-ray and optical 
observations (mass transfer rate, position on the CM dia- 
gram, characteristic ages of the parent stellar cluster and 
the surrounding bubble nebula), the system would be con- 
sistent with a ~ 50 — IOOMq black hole accreting from a 
12-15M© star that fills its Roche lobe at terminal age main 
sequence. At this stage of the binary evolution, such a star 
would have a radius of ~ 8 — 9 x 10 11 cm and a separa- 
tion from the BH of ~ 4 x 10 12 cm, and so it will intercept a 
fraction ~ 10— 15% of the X-ray photons assuming isotropic 
emission. Note, however, that the fraction of photons inter- 
cepted by the disc exceeds that from the donor star (which 
is also partly screened by it). 

Finally, the detection of the orbital modulation would 
open the way to perform a direct dynamical measurement 
of the BH mass in NGC 1313 X-2, the first time ever for 
a ULX, although a recent Gemini sp ectroscopy campaig n 
failed to reveal regular modulations (jRoberts et al.ll2010l ). 
The binned light curve suggests that the periodicity may be 
there, but the low statistical significance of the B band mod- 
ulation, along with the lack of detection in the V band, make 
its identification uncertain. A dedicated photometric mon- 
itoring campaign under homogeneous observing conditions 
to minimize systematic uncertainties are needed to confirm 
it. 



ACKNOWLEDGMENTS 

DI and LZ acknowledge financial support from INAF 
through grant PRIN-2007-26. RS acknowledges financial 
support from a Curtin University senior research fellowship, 
and hospitality at the Mullard Space Science Laboratory 
(UCL) and the University of Sydney during part of this 
work. We thank the referee for his/her very careful reading 
of the manuscript and the valuable comments that helped 
us to improve significantly upon the previous version of our 
paper. This paper is based on observations collected at ESO 
VLT. Observations were also used from NASA/ESA HST, 
obtained from the data archive at the Space Telescope In- 
stitute. STScI is operated by the association of Universities 
for Research in Astronomy, Inc., under the NASA contract 
NAS 5-26555. 



REFERENCES 

Bochkarev, N. G., Karitskaia E. A., Shakura N. I., 1979, 

Astron. Zh., 56, 16 
Colbert, E. J. M., & Mushotzky, R. F., 1999, ApJ, 519, 89 
Colbert, E. J. M., & Ptak, A. F., 2002, ApJS, 143, 25 
Copperwheat, C, Cropper, M., Soria, R., & Wu, K. 2007, 

MNRAS, 376, 1407 



8 L. Zampieri, D. Impiombato, R. Falomo, F. Grise and R. Soria 



Davies, S. R., 1990, MNRAS, 244, 93 

Dolphin, A.E., 2009, PASP, 121, 655 

Fabbiano, G. 2006, ARA&A, 44, 323 

Feng, H., Kaaret, P., 2006, ApJ, 650, 75 

Feng, H., Kaaret, P., 2008, ApJ, 675, 1067 

Feng, H., Soria, R. 2011, New Astronomy Review, in press 

Grise, F., et al., 2006, IAU Symp. 230, p. 293 

Grise, F., et al., 2008 , A&A, 486, 151 

Grise, F., et al., 2009, AIPC, 1126, 201 

Kaaret, P., Ward, M. J., Zezas, A. 2004, MNRAS, 351, L83 

King, A. R., Davies, M. B., Ward, M. J., Fabbiano, G., & 

Elvis, M., 2001, ApJ, 552, 109 
King, A. R., 2008, MNRAS, 385, L113 
Liu, J.-F., & Bregman, J. N., 2005, ApJS, 157, 59 
Liu, Q. Z., & Mirabel, I. F. 2005, A&A, 429, 1125 
Liu, J. et al., 2004, ApJ, 602, 249 
Liu, J., et al., 2007, ApJ, 661, 165 
Liu, J., et al., 2009, ApJL, 690, L39 

Mapelli, M., Colpi, M., & Zampieri, L. 2009, MNRAS, 395, 
L71 

Mendez, B., et al., 2002, AJ, 124, 213 
Mucciarelli, P., et al, 2005 , ApJL, 633, L101 
Mucciarelli, P., et al., 2007, ApJ, 658, 999 
Pakull M. W., Mirioni L., 2002, New Visions of the X-ray 
Universe in the XMM-Newton and Chandra Era, 26-30 
November 2001, ESTEC, Noordwijk (0202488) 
Pakull, M.W., et al., 2006, IAU Symp., 230, 293 
Patruno, A., Zampieri, L., 2008, MNRAS, 386, 543 
Patruno, A., Zampieri, L., 2010, MNRAS, 403, L69 
Ramsey, C.J., et al., 2006, ApJ, 641, 241 
Roberts, T. P., & Warwick, R. S., 2000, MNRAS, 315, 98 
Roberts, T. P., Levan, A. J., Goad, M. R., 2008, MNRAS, 
387, 73 

Roberts, T. P., Gladstone, J. C, Goulding, A. D., Swin- 
bank, A. M., Ward, M. J., Goad, M. R., Levan, A. J., 

2010, AN, in press (|arXiv:1011.2155|l 
Rizzi, L., et al., 2007, ApJ, 661, 815 
Saesen, S., et al., 2010, A&A, 515, id.A16 
Soria, R., et al., 2005, MNRAS, 356, 12 

Swartz, D. A., Ghosh, K. K., Tennant, A. F., & Wu, K., 
2004, ApJS, 154, 519 

Swartz, D.A., Soria, R., Tennant, A.F., & Yukita, M., 2011, 
ApJ, in press l|arXiv:1108.1372|l 

Tully, R. B., 1988, Nearby Galaxies Catalog. Cambridge 
Univ. Press, Cambridge 

Uslenghi, M., Falomo, R., 2008, in Modelling and Simula- 
tion in Science, ed. V. Di Gesu et al. (Hackensack, NJ: 
World Scientific), 313 

van Paradijs, J., McClintock, J. E. 1995, in X-ray binaries, 
Edited by W. H. G. Lewin, J. van Paradjis, E. P. J. van 
den Heuvel, Cambridge: Cambridge University Press, p. 
58 

Walton, D. J., Roberts, T. P., Mateos, S., & Heard, V. 

2011, MNRAS, 1147 

Zampieri, L., et al., 2004, ApJ, 603, 523 

Zampieri, L., & Roberts, T. P. 2009, MNRAS, 400, 677 



